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I.  INTRODUCTION 


1  la 

A  previous  paper  *  examined  the  harmonic- oscillator  model  of 
2  3 

Drickamer  et  al *  for  the  effects  of  high  pressure  on  molecular,  elec- 

4 

tronic  spectra.  The  parameterization  of  the  model  by  Okamoto  et  al 

was  reformulated  to  conform  with  the  conclusions  of  Lin5  and  tested  on 

the  data  of  Okamoto  and  Drickamer.6  Although  the  reformulation  im¬ 
proved  the  agreement,  the  results  were  not  entirely  satisfying.  As  one 
of  the  possible  causes  for  discrepancies,  I  suggested  anharmonicity. 

For  the  next  level  of  sophistication  in  this  direction,  the  incorpora¬ 
tion  of  the  Horse  potential  suggests  itself. 

The  present  paper  derives  expressions  analogous  to  those  of  Refer¬ 
ence  1,  based  upon  the  Morse  potential  rather  than  the  harmonic  oscil¬ 
lator.  Since  Lin5  has  shown  that  his  quantum-mechanical  treatment  leads 
to  the  Drickamer  equations  when  the  same  assumptions  are  made,  we  shall 
follow  Drickamer 's  simpler  and  more  transparent  approach  here.  Part  II 
summarizes  the  approach  and  introduces  the  basic  equations  that  will  be 
required.  The  expressions  for  peak  height  are  derived  in  Part  III,  and 
those  for  bandwidth  in  Part  IV.  Part  V  casts  the  expressions  in  a 
compact  form  suitable  for  experimental  testing,  and  Part  VI  makes  con¬ 
cluding  remarks.  The  Appendix  lists  several  expansions  in  Maclaurin 
series  that  are  used  in  the  text. 


II.  GENERAL  CONSIDERATIONS 

We  assume  that  pressure  couples  to  a  Morse  potential  along  a  single 
configuration  coordinate,  viz.,  volume.  The  potential  functions  of 
the  ground  and  excited  electronic  states  are  then  given,  respectively, 

by  2 

V  =  D(1  -  e‘SQ)  +  cpQ  (2.1) 

atld  V*  =  d*(1  -  e'6*^)2  +  Tcpq  *  E0*  (2*2) 

TI  Tompkins,  R.C.,  J.  Chem.  Phys.  69,  579  (1978). 
la.  Tompkins,  R.C.,  Effects  of  High  Pressure  on  Molecular  Electronic 
Spectra;  A  Comparison  of  Models,  BRL  R  No.  2014,  (1977) . (A045414) 

2.  Drickamer,  H.G.,  Frank,  C.W.,  and  Slichter,  C.P.,  Proc.  Nat.  Acad. 
Sci.  (USA)  69,  933  (1972). 

3.  Drickamer,  H.G.,  and  Frank,  C.W.,  "Electronic  Transitions  and  the 
High  Pressure  Chemistry  and  Physics  of  Solids,  "London:  Chapman 
and  Hall  (1943). 

4.  Okamoto,  B.Y.,  Drotning,  W.E.,  and  Drickamer,  H.G.,  Proc.  Nat. 

Acad.  Sci.  (USA)  71,  2671  (1973). 

5.  Lin,  S.H.  J.  Chem.  Phys,  59,  4458  (1973) 

6.  Okamoto,  B.Y.,  and  Drickamer,  H.G.,  J.  Chem.  Phys.  61,  2870  (1974). 


where  D,  0,  D*,  e*,  c  are  constants,  p  is  pressure,  Q  is  the  configura¬ 
tion  coordinate,  r  is  the  ratio  of  the  coupling  constant  in  the  excited 
state  to  that  in  the  ground  state,  and  Eq  and  q  are  the  energy  and  co¬ 
ordinate  separations,  respectively,  between  the  potential  minima.  Note 
that  when  the  Morse  function  is  reduced  to  the  harmonic  function,  we 

2  2 
obtain  a  of  References  1-6  as  20  D. 

We  shall  obtain  an  expression  for  the  optical  absorption  energy  by 
projecting  the  minimum  of  eq.  (2.1)  upon  eq.  (2.2)  in  accordance  with 
the  Franck-Condon  principle,  and  mutatis  mutandis  for  emission.  We  shall 
obtain  the  bandwidths  by  projecting  the  terminal  points  of  oscillation 
of  the  lowest  vibrational  level. 

Differentiation  of  eqs.  (2.1)  and  (2.2)  yields 

3V/3Q  =  20De”8Q(l  -  e"BQ)  +  cp  (2.3) 

and 

3V*/3Q  =  2Rp0Dye‘p6Q(l  -  Ye'pBQ)  +  rep,  (2.4) 

where  R  ■  D*/D,  p  ■  0*/0,  and  y  ■  epB<*  •  Setting  these  expressions 
equal  to  zero  and  solving  as  quadratic  equations  in  the  exponentials, 
we  obtain 


e~B<*m  *  h(l  ♦  /l+2cp/0b  )  (2.5) 

and 

e'pBQ*m  «  (1  +/  l+2rcp/Rp0D)/2y.  (2.6) 

Hence 

On  =  [log  2  -  log(l  +  /l+2cp/0D) ]/3  (2.7) 

and 


Q£  =  q  +  [log  2  -  log(l  +  A+2rcp/Rp0D)]/p0.  (2.8) 


The  potential  minima  are  then 


v.  =  wd .  /i72cp78p) *  1  ♦  m»s  y.  i»i»  *  (2.9) 

V*  =  *sRD(l  -  A+2rcp/Rp0D)  +  crqp  +  E 

°  (2.10) 

+  1  *  2 [log  2  -  log(l  ♦  2rcp/R  D)] 

2p0  rcp* 
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It  will  be  convenient  to  obtain  at  this  point  certain  other  quan¬ 
tities  that  will  be  required  in  the  sequel.  Differentiation  of  eqs. 
(2.9)  and  (2.10)  with  respect  to  pressure  yields 


V^(P)  »  2fjDQjje_8Qm(l  -  e”8^)  +  c^p  ♦  QJ 

(2.11) 

and 

(p)  =  2RpBDYQ|J/e'paQm(l  -  Ye"p8Q*m)  +  rc  (QjJ' JwQ^) . 

(2.12) 

Evaluating  eqs.  (2.7)  -  (2.12)  at  p»0,  we  obtain 

O 

w 

u 

o 

(2.13) 

•< 

o 

v_-/ 

II 

O 

*» 

(2.14) 

to)  -  o. 

(2.15) 

/ — \ 

o 

w 

n 

(2.16) 

Vm^  =  Eo» 

(2.17) 

and 

Vm  "  rcq* 

(2.18) 

(We  shall  not  require  Q^(0)  or  Q£'  (0).) 


III.  PEAK  HEIGHT 

A.  Absorption. 

We  project  the  potential  sdnimum  of  the  ground  state  upon  the 
potential  function  of  the  excited  state  by  substituting  eq.  (2.7)  into 
eq.  (2.2): 

V*(V  *  RD^  “  2y(1  +  v'i+2cp7eD)P/2P+  y2(l  +  /l+2cp/BD)2p/22p] 

(3.1) 

+  rc[log  2  -  log(l  +  /l+icp/BD)]p/8  +  EQ. 


Subtraction  of  eq.  (2.9)  yields  the  optical  absorption  energy. 


It  will  be  convenient  to  obtain  at  this  point  certain  other  quan¬ 
tities  that  will  be  required  in  the  sequel.  Differentiation  of  eqs. 
(2.9)  and  (2.10)  with  respect  to  pressure  yields 


V;(P)  -  20DQjie'e<lm(l  -  e_$Qm)  +  0(0^  ♦  QJ 

(2.11) 

and 

(p)  =  2Rp$DyQ|J7e'peQm(l  -  ye"p$Q*m)  +  rc (Q*  f»Q£) . 

(2.12) 

Evaluating  eqs.  (2.7)  -  (2.12)  at  p»0,  we  obtain 

(^(0)  «  0, 

(2.13) 

O 

II 

r~\ 

O 

(2.14) 

^(°)  a 

(2.15) 

0^(0)  -  q. 

(2.16) 

V*(0)  =  B0. 

(2.17) 

and 

V*7 (0)  =  rcq. 

(2.18) 

(We  shall  not  require  Q^(0)  or  Q£7  (0).) 

III.  PEAK  HEIGHT 

A.  Absorption. 

We  project  the  potential  minimum  of  the  ground  state  upon  the 
potential  function  of  the  excited  state  by  substituting  eq.  (2.7)  into 
eq.  (2.2): 

V*(Qm)  =  RD[1  -  2y(l  +  /l+2cp/0D)p/2p+  y2(l  +  /l+2cp/0D) 2p/22p] 

(3.1) 

+  rc[log  2  -  log(l  +  /l+2cp/0D)]p/0  +  Eq. 

Subtraction  of  eq.  (2.9)  yields  the  optical  absorption  energy. 
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ea  “  V*<V  -  Vm 

-  RDy2-p (1  +  ✓l+2cp/gD)p[2_PY(l  +  /l+2cp/6D)p  -  2]  +  K 

(3.2) 

+%D/l+2cp/3D  -  cp/20  +  (cp/0) (1-r) [log(l  +  /l+2cp/gD)  -  log  2], 


where  K  3  E0  +  D(R  -  H) t  which  is  independent  of  pressure. 


The  experiments 
ing  evidence  that  6 

Si 


of  Okamoto  and  Drickamer^*^  have  provided  convinc- 
is  a  quadratic  function  of  p.  Therefore,  we  shall 


expand  eq.  (3.2)  in  a  Maclaurin  series  to  second  order  in  p,  using 
eqs.  (A. 1),  (A. 2),  and  (A. 3).  After  subtracting  the  energy  at  zero 
pressure,  we  obtain  the  pressure  shift  of  the  absorption  peak. 


6  t 

a 


Rpx(x.-A)cp  + 

P 


Rpy[p (2y  -  1)  -  3(y  -  2)]  +  (1  -  2r)  2  2 

T-  C  P 

43  D 


(3.3) 


B.  Emission. 

We  project  the  potential  minimum  of  the  excited  state  upon  the 
potential  function  of  the  ground  state  by  substituting  eq.  (2.8)  into 
eq.  (2.1): 

V(Q*)  «  D[1  -  2(*  *  l/l+2rcp/RpeD)  1/p  +  fl  +  /l+2rcp/RpBP)  2/p^ 

tB  *y  2y 

(3.4) 

+  cpq  +  (cp/p0)  [log  2  -  logfl  +  /l+2rcp/RpBD) ]  . 

Subtraction  from  eq.  (2.10)  yields  the  optical  emission  energy, 
g  =  ^-1/p^l  ♦  ✓H-2rcp/RpBD)l/p^2  _  ^1  +  /l+2rcp/Rp0Dj  1/p ^ 

-  *iRD/l+2rcp/Rp0D  +  rcp/20D  +  (r  -  l)cqp  (3.5) 

+  (cp/pB) (r  -  1)  [log  2  -  log (1  +  /l+2rcp/Rp8D)]  +  K7  , 

T.  Okamoto,  B.Y.,  and  Drickamer,  H.G.,  J.  Chem.  Phys.,  61,  2870  (1974). 
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where  K7  ■  EQ  +  (R/2  -  1)D,  which  is  independent  of  pressure.  Again  . 

we  expand  in  Maclaurin  series  to  second  order  in  p,  using  eqs.  (A. 4), 
(A. 5),  and  (A. 6),  and  subtract  the  energy  at  zero  pressure.  The  pres¬ 
sure  shift  of  the  emission  peak  is  then. 


St  -  n'1/ll(i  -  y'1/o)  «  (r  -  l)  «P2eqcp 
e  Rp2s 


(3.6) 


.  r  Y~1/p{r[(l  -  2y  ~1/p  -  3p(l  -  y'1/p)1  ♦  R(2-r)>  2  2 

p  . 

4R  p^D 


IV.  BANDWIDTH 


A.  Absorption. 


Let  us  represent  the  terminal  points  of  oscillation  in  the  lowest 
vibrational  state  of  the  electronic  ground  state  by  Q+  and  Q  (collec¬ 
tively,  .  Then  the  bandwidth  is  given  by  the  energy  difference 

between  their  projections  upon  the  potential  function  of  the  excited 
state: 


E*  -  E*  -  RD[  (1-Y<)2  -  (1-Y^)2]  ♦  rc(Q+-Q_)p, 

where 

-3Q+ 

♦±  *  e 


(4.1) 


(4.2) 


It  will  be  convenient  to  rearrange  eq.  (4.1)  in  the  form 


E*-E*  «  RDy[y(42P-*?P)  -  2(*J  -  4P)1  ♦  rc(Q+-Q_)p.  (4.3) 

Following  the  experimental  evidence  of  References  6  and  7  that  the  band 
width  is  a  linear  function  of  pressure,  let  us  expand  the  bandwidth 
in  a  Maxlaurin  series  to  first  order  in  p: 


E*-E*  =  (E*-E*)0  +  (E*-E*)^p, 


(4.4) 


where  a  zero  subscript  denotes  evaluation  at  p»0,  and  a  prime  denotes 
a  derivative  with  respect  to  p.  Differentiation  of  eq.  (4.3)  yields 


/ 


(E*  .  E*)'  «  2RD»>y[y(*JP"V+  -  *?P‘V_)  -  C*J'V+  -  *P'Vj] 

(4.S) 


+  rc(Q+  -  Q_)  +  rc(Q^  -  Q^ip. 


The  energy  at  the  terminal  points  of  oscillation  cai  be  expressed 

by  2 

D(1  -  *±r  +  CQ±P  =  Vm  +  kT/2,  (4.6) 

where  k  is  the  Boltzmann  constant  and  T  is  temperature.  Although  this 
transcendental  equation  cannot  be  solved  analytically  for  Q+,  we  can 

evaluate  it  at  p=0  with  the  aid  of  eq.  (2.14).  Defining  a=  /kT/2D, 
we  have 


1  -  <fr±(0)  =  ±o 

(4.7) 

or 

♦+(0)  a  1  *  « 

(4.8) 

and 

Q±(0)  =  -log(l  ?o)/B  . 

(4.9) 

The 

signs  have  been  chosen  such  that  Q+(0)>0  and  Q_(0)<0. 

to  p 

The  implicit  differentiation  of  eqs.  (4.2)  and  (4.6) 
yields 

,  .  -W± 

♦i  * 

with  respect 

(4.10) 

and 

-2DV±(1  -  $±)  +  cQ'±p  +  cQ±  =  Vj  . 

(4.11) 

Recalling  eq.  (2.15),  we  obtain 

v  (0)  =  +  c  W1  L.  ZL 

J  2a£5D 

(4.12) 

and 

.  <(°) 

Qi(0)  =  -  — —  . 

e<^+(0) 

(4.13) 
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Evaluation  of  eq.  (4.3)  and  (4.5)  at  p=0  introduces  complicated 
functions  of  a.  Since  the  energy  at  the  potential  minimum  is  certainly 
much  less  than  the  dissociation  energy,  D,  a  «  1.  Therefore,  we  can 
simplify  by  means  of  Maclaurin  series  in  a  (eqs.  (A. 7),  (A. 8),  and 
(A. 9)),  retaining  only  the  first  nonvanishing  power  of  o.  Thus,  we 
obtain 


(E*  -  E_)0  =  4RDayp (1  -  y)  .  (4.14) 

and 

(E*  -  E*)'0  =  (Rcayp/8)  [2p  -  3  -  y(4p  -  3)]  +  2rca/0.  (4.15) 

Substituting  eqs.  (4.14)  and  (4.15)  into  eq.  4.4,  we  have 

(4.16) 

E*  -  E*  a  4RDayp (1  -  y)  +  (ca/0)  {Ryp[2p-3  -  y(4p-3)]  +  2r}p. 

The  Gaussian  half-width  is  obtained  by  multiplying  by  (2  log2),*: 

e  =  4Ryp  (1-y)  (kTD  ^2)**  +  [c(kT  log2)Ss/0/D]  {Ryp  [2p-3-y  (4p-3)  ] 

a 


+  2r}p. 

B.  Emission. 


(4.17) 


Analogous  to  eq.  (4.1)  we  have 

E+  -  E_  =  D[(l  -  **1/p)2  -  (1  -  *_1/p)2]  +  c(Q*  -  Q*)p  , 

where 

-p0q: 

K  -  e 


(4.18) 

(4.19) 


Rewrite  eq.  (4.18)  as 

E+  -  E_  =  D[(4*2/p  -  **2/p)  -  2(4*1/P  -  4!1/P)]  ♦  c(Q*  -  Q*)p. 
Then 

(E+  -  E _)'  =  2 Dp1  ((♦*2/p-14j/  -  *!2/p’14!/)  -  (♦:1/P_V 
-  ♦;1/p_1^)]  +  c(Q*  -  Q*)  ♦  c (Q*;  -  Q*')p. 


(4.20) 


(4.21) 
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The  energy  at  the  terminal  points  of  oscillation  in  the  excited 
is 

RD(1  -Y*i)2  ♦  rcQ*p  +  Eq  -  VJ;  ♦  kT/2. 

Defining  a*  =  /kT/2RD  and  recalling  eq.  (2. 17),  we  have 

1  -y$|(0)  =  ±  a* 
or 

+i(0)  =  (1  r  a*)/Y 

and 

QJ(0)  -  q  -  lpt(p-e-  ■a-*3-  • 


Furthermore, 

-peo: 

K'  =  -P&Q+'  e 

and 


-2RDy4>^/  (1  -Y ♦*)  ♦  re  (Q*7  p  ♦  Q*)  =  V*'  , 


whence  with  eq.  (2.18) 


and 


QJ'  (0)  =  I  - 

2RpVd«*(1  *  a *) 


Applying  Maclaurin  series  as  before,  we  have 


E+  -  E 


.  -  V+  -  EJ0  +  <E+  ~  E-)'op 

«.  !l°  +  2}p 

py2/p  pB  Rpy27^ 


state 

(4.22) 

(4.23) 

(4.24) 

(4.25) 

(4.26) 

(4.27) 

(4.28) 

(4.29) 

(4.30) 
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J 


and  the  Gaussian  half-width, 


e  *  r(kTD  log  2)ls(Y1/el)/PY2/(VR 

G 


+ 


c_  AT  lot  2Atr[(»-2p)T1/p- 

‘  RD  J 


+  2>p  . 


V.  PANSYNTHESIS 


(4.31) 


It  remains  to  cast  the  model  in  a  form  convenient  for  application 
to  experimental  data.  Following  the  procedure  adopted  in  Reference  1, 
we  square  eqs.  (4.17)  and  (4.31)  and  subtract  the  terms  independent  of 
pressure  to  obtain 


6ef/N2  =  (c/0)Ryp(1-y)  (Ryp  [2p-3-Y(4p-3) ]  +  2r>p 

ci 


+  (c2/802D)(RYP[2p-3-Y(4p-3)]  +  2r}2p2 


(5.1) 


and 


6ef/N2  *  (c/Rp20Y2/P)(r((2-2p)Y1/p-(4-3p)]/RpY2/p  +  2}p 


+  (c2/8RDp202)(r[(2-2p)Y1/p-(4-3p)]/RPY2/p  +  2}2p2, 
where,  as  in  Ref.  1,  N  =  (8kT  log  2)^. 


(5.2) 


Spectroscopic  measurements  such  as  those  of  References  6  and  7 
can  be  expressed  in  the  form 


and 


66a  *  alp  +  a2p2,  (5-3) 
6ee  s  elp  +  e2p2,  (5-4) 
6ea  a  dap’  (5’5) 

6ee  =  dep>  (5-6) 


where  the  coefficients  on  the  right  are  empirical  fitting  parameters. 


IS 


As  was  pointed  out  in  Ref.  1,  we  can  write 

fie^/N2  *  2Eap  +  Dap2  (5.7) 

and 

6e2/N2  =  2Egp  +  Dep2  ,  (5.8) 

where  E  =  d  e°/N2,  D  *  d2/N2,  e°  is  the  absorption  half-width  at 

Ado  a  a  a 

zero  pressure,  etc. 

Equations  (5.1)  -  (5.8)  with  eqs.  (3.3)  and  (3.6)  give  us  finally, 
defining  8  *  8/c, 

*2  =  RPY(Y  -  D/B  >  (5.9) 

a 2  «  (Rpy[p(2Y-1)  -  3(y-2)]  +  (l-2r) }/482D,  (5.10) 


Ea  =  Rpy(1-y)(Rpy[2p-3-y(4p-3)]  +  2r)/2$,  (5.11) 

Da  =  {Rpy[2p-3-y(4P-  3)]  +  2r>2/882D,  (5.12) 

ej  =  rY"1/p(l-Y‘1/p)/Rp2B  ♦  q(r-l),  (5.13) 

e2  =  rY"1/p{r[(l-2Y"1/p-3p(l-Y"1/p)]  ♦  R(2-r) >/4R2p2  g2D,  (5.14) 
Ee  -  (r[(2-2p)Y1/5(4-3p)]/RpY2/p  +  2}/2Rp2$Y2/?  (5.15) 

De  =  (r[(2-2p)Y1/p-(4-3p)]/RpY2/p  +  2)2/8RDp2g2  .  (5.16) 


VI .  CONCLUDING  REMARKS 

Proceeding  in  the  spirit  of  Drickamer  et  sU  ',s'  ,  we  have  derived 
from  a  Morse  potential  a  set  of  equations  for  evaluating  the  configura¬ 
tion-coordinate  parameters  that  describe  the  effects  of  high  pressure 
on  optical  absorption  and  emission  spectra.  These  eight  equations 
relate  seven  unknown  parameters  (6,  D,  q,  R,  r,p)  to  eight  experimentally 
accessible  quantities.  We  have  made  the  following  assumptions: 
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1.  The  pressure  dependence  of  peak  height  is  second  order,  and 
that  of  bandwidth  first  order. 

2.  The  constants  in  the  Morse  function,  6  and  D,  are  independent 
of  pressure. 

3.  The  parameters  can  be  regarded  as  averages  over  all  normal 
modes  and  configuration  coordinates. 

The  first  assumption  is'  consistent  with  experimental  evidence**  ’ ^ .  The 

second  assumption  is  suggested  by  Lin's5  conclusions  from  first-order 
perturbation  theory  in  the  harmonic  approximation.  The  third  assump¬ 
tion  is  necessary  to  avoid  an  intractable  proliferation  of  parameters. 

If  future  experimental  results  should  force  a  relaxation  of  the  first 
assumption,  it  is  necessary  only  to  carry  the  Maclaurin  expansions  to 
higher  order. 

Resources  are  not  at  present  available  to  exercise  eqs.  (5.9)  - 
(5.16)  on  the  experimental  data.  It  might  even  be  questioned  whether 
such  an  endeavor  is  worthwhile  until  more  data  are  at  hand.  In  prin¬ 
ciple  the  Morse  potential  should  be  more  accurate  than  the  harmonic 
oscillator.  It  would  be  desirable  to  have  a  perturbation  treatment 

similar  to  that  applied  by  Lin5  for  the  harmonic  oscillator. 

8  9 

The  significant  work  by  Berry  and  Williams  and  Curie  and  Williams 
in  deriving  a  model  of  pressure  effects  from  first  principles  has  come 
to  my  attention  since  completing  the  work  presented  here.  Although 
their  model  in  its  harmonic  approximation  works  well  for  inorganic 
crystals,  it  is,  like  the  original  Drickamer  model,  qualitatively  incon¬ 
sistent  with  experiments  on  aromatics.  However,  a  Morse  treatment  of 
that  model  along  the  lines  of  the  present  paper  may  bear  fruit. 

Reference  la  details  the  motivation  of  this  work  from  an  interest 
in  pressure  effects  on  burning  rates  of  propellants.  It  has  now  be¬ 
come  clear  that  the  effects  under  study  are  unlikely  to  change  dissocia¬ 
tion  energy  sufficiently  to  influence  reaction  rates.  Since  this 
approach  to  propellant  phenomenology  no  longer  appears  promising,  the 
work  has  been  discontinued. 


8.  Berry,  D.,  and  Williams,  F.,  Acta  Physica  Polomica  A54,  705  (1978). 

9.  Curie,  D.  and  Williams,  F.,  in  Luminescence  of  Inorganic  Solids, 

B.  Di  Bartolo  ed.,  New  York:  Plenum  Press  (1978). 
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This  appendix  gives  the  requisite  expansions  in  Maclaurin  series. 


/l+2cp/0D  =  1  +  cp/fJD  -  c2p2/2eV  ♦  ...  (A. I) 
log(l  +  ^l+2cp/0D)  =  log  2  ♦  cp/20D  -  3c2p2/402D2  +  ...  (A. 2) 
[(1  +  /l+2cp/SD)/2]p  =  a  +  pcp/20D  +  p (p-3)c2p2/802D2  +  ...  CA.3) 
/l+2rcp/Rp0D  =  1  +  rcp/RpgD  -  r2c2p2/2R2p202D2  ■*■...  (A. 4) 


log(l  +  »/l+2rcp/Rp0D)  =  log  2  +  rcp/2Rp0D  -  3r2c2p2/4R2p202D2  +  ...  (A. 5) 

[(1  +  *'l+2rcp/Rp0D)/2] 1/p  =  1  +  rcp/2Rp2$D 

(A.uj 

♦  r2c2(l-3p)p2/8R2p4B2D2  +  ... 

log (1  ±  a)  =  ±  a  -  Jja2  ±(l/3)a3  -  ...  (A. 7) 

(1  ±  a)CT  =  1  ±  aa  +  !sa(ar-l)a2  ±[a(a-l)  (<J-2)]a3/6  +  ...  (A. 8) 

(1  ±  a)alog(l  ±a)  =  ±  a  +  H(.2o-l)a2  ±  [3a(o-2)  +  2]a  3/6  +...  (A. 9) 
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